Quantum Hall Effect in Thin Films of Three-Dimensional Topological Insulators 
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We show that a thin film of a three-dimensional topological insulator (3DTI) with an exchange field 
is a realization of the famous Haldane model for quantum Hall effect (QHE) without Landau levels. 
The exchange field plays the role of staggered fluxes on the honeycomb lattice, and the hybridization 
gap of the surface states is equivalent to alternating on-site energies on the AB sublattices. A peculiar 
phase diagram for the QHE is predicted in 3DTI thin films under an applied magnetic field, which 
is quite different from that either in traditional QHE systems or in graphene. 
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Introduction. -Haldane was the first to realize that 
quantum Hall effect (QHE) can happen in systems with- 
out Landau levels (LLs) . He proposed a spinless electron 
model defined on a two-dimensional (2D) honeycomb 
lattice with staggered magnetic fluxes and alternating 
on-site energies on the AB sublattices. In the model, the 
electrons retain their usual Bloch state character, but ex- 
hibit a nonzero quantized Hall conductivity Uxy by simu- 
lation of the "parity anomaly" of (2 + 1)D field theories, 
now termed as quantum anomalous Hall effect. The Hal- 
dane model breaks time-reversal symmetry. Kane and 
Mele suggested that the single-atomic-layer graphene 
with intrinsic spin-orbit coupling constitutes two copies 
of the Haldane model. Since the two spin components are 
mutually conjugate under time reversal, the Kane-Mele 
model restores the time reversal symmetry. While the 
two spin components contribute oppositely to the charge 
Hall conductivity, the model displays a novel quantum 
spin Hall effect (QSHE). The QSHE has been experimen- 
tally observed in quantum wells of mercury telluride Q , 
instead of graphene. A great deal of research interest in 
the QSHE in recent years has led to the birth of topolog- 
ical insulators, which include a class of 2D and 3D band 
insulators with nontrivial band topology 

A 3D topological insulator 0, i (3DTI) is featured 
with the presence of conducting surface states in the bulk 
band gap, which are topologically protected and insen- 
sitive to impurity scattering. On an isolated surface of 
the 3DTI, the energy spectrum of the surface states con- 
stitutes a single massless spin-momentum-locked Dirac 
cone I?], which offers a unique platform to realize some 
exotic relativistic quantum phenomena, such as Majo- 
rana fermions [8]. In a perpendicular magnetic field Bq, 
the surface states are quantized into LLs with energies 
En — sgn(n) yj2eBQv'^\n\ (in fi — c — 1 units), where 
n = 0, ±1, ±2 • • • and vp is the Fermi velocity. Due to the 
existence of the zero-mode (n — 0) LL, which is particle- 
hole symmetric, the Hall conductivity is expected to be 
half-integer quantized Q axy — Vj^, where v — i \ 
with i being an integer. Exciting experimental progress 
has been made recently. Cheng et al. and Hanaguri et 



al. |9| detected the surface states in Bi2Se3 to form LLs 
in the presence of an applied magnetic field. Very re- 
cently, Briine et al. [lo| carried out a Hall conductivity 
measurement in a high-quality strained HgTe layer, and 
observed a sequence of Dirac-like Hall plateaus, which 
was ascribed to the topological surface states. It was 
also suggested theoretically that electrons in the surface 
states of a 3DTI with exchange splitting could give rise 
to a QHE in the absence of an applied magnetic field. 



i.e., the quantum anomalous Hall effect For a thin 
film of a 3DTI, the finite-size confinement will mix the 
surface states at the top and bottom surfaces and create 
a hybridization energy gap in the spectrum of the surface 
states I12I . |l3| . The hybridization gap can be controlled 



by tuning the thickness of the 3DTI film [H [H. The 
exchange field (or equivalently a Zeeman energy) and the 
hybridization gap effectively add different types of mass 
to the Dirac fermions [13] and provide ways to control 
the properties of the surface states, but their effect on 
the QHE has not been well understood. 

The Haldane's QHE model '1'] has been difficult to be 
realized in condensed matter systems, though Shao et 
al. suggested recently that it could be simulated by 
using ultracold atoms. In this Letter, we propose that a 
thin film of a 3DTI in the presence of an exchange field is 
a natural realization of the Haldane model. The exchange 
field plays the role of the staggered fluxes in the Haldane 
model, and the hybridization gap is equivalent to the 
alternating on-site energies on AB sublattices. Under an 
applied magnetic field, a 3DTI thin film and the Haldane 
model are still equivalent in Hamiltonian form, so that 
we can study the QHE in the 3DTI film conveniently 
based upon the Haldane model. The other purpose of 
this Letter is to report a peculiar phase diagram for the 
QHE in thin films of 3DTIs, due to the interplay between 
Zeeman energy g and hybridization gap A. For A = 0, 
a nonzero g only shifts the positions of the LLs without 
creating new Hall plateaus. If 17 = but A ^ 0, a new 
1^ = plateau will appear, in addition to the original v = 
2(£+i) odd-integer plateaus. The simultaneous presence 
of nonzero g and A causes splitting of the odd-integer 
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Hall plateaus, and so all integer = £) plateaus emerge. 
More interestingly, as the product of g and A is equal 
to certain critical values, the split plateaus can merge 
again, and most odd-integer (or even-integer) plateaus 
disappear. 

Thin Films of SDTIs.-We start with an effective 
Hamiltonian for electrons in the surface states of a 3DTI 
thin film [1^ with an exchange energy 11 1 



H (k) = -Dk^+VF {ky&x - kxay) + { — -Bk'^)Tzaz+gffz- 

(1) 

Here, the first two terms are the Hamiltonian for two 
isolate surfaces with vf the Fermi velocity, k is the 
momentum with respect to Dirac point F, and iTq are 
the Pauli matrices for electron spin with a = x,y,z. 
The third term stands for the coupling between the 
upper and lower surface states of the thin film with 
Tz = 1 (—1) representing the bonding (antibonding) be- 
tween them, and the last one is the exchange field with 
strength g. We will consider the case, where the con- 
dition — > for band inversion is satisfied [l3| . 
Since Hamiltonian (1) is block diagonal for = ±1, the 
Chern numbers Cr^ of the occupied valence bands with 
Tz = +1 and —1 can be calculated separately, yielding 
Cr, = ~Tz [sgn(B) + sgn(A + 2gTz)] /2. The total Chern 
number C = C+ + C_ is given by 
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C=-[sgn(A-25)-sgn(A + 25)] 



(2) 



We note that even though the r^-dependent Chern num- 
bers depend on the sign of parameter i?, the total Chern 
number is independent of i?. In Fig. 1, we show the phase 
diagram of the 3DTI thin film in the presence of exchange 
splitting. Lines g = —A/2 and g = —A/2 are the criti- 
cal boundaries between the normal insulator with (7 = 
and the Chern insulator with C = ±1, the latter being 
in the lower and upper shaded regions. In the regions of 
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C = ±1, quantum anomalous Hall effect with axy = ±x 
occurs. 
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FIG. 1: Schematic phase diagram on the hybridization po- 
tential A and exchange energy g plane, in which two solid 
lines g = ±A/2 separate four phase regions with C = ±1 and 
C = 0. 



Haldane Model.-The Hamiltonian of the Haldane 



model reads Q 

H^-J2 tele, + J2 ^26^*"^ cjc, + J2 y^4c^ ■ (3) 

Here, c\ (q) is the fcrmion creation (annihilation) op- 
erator at site i, t [12) is the hopping integral between 
the nearest-neighbor (next-nearest-neighbor) sites i and 
j, and Vi = ±M for sublattices A and B, respectively. 
(/>y = ±(/) is the accumulated Peierls phase from site i 
to its second neighbor j, due to a staggered magnetic- 
fiux density, where the positive sign is taken for an elec- 
tron hopping along the arrows indicated in the Fig. 1 
of Ref. [1]. It has been shown fl]] that the Haldane 
model (3) exhibits a quantized Hall conductivity ±^ for 
3-\/3|t2 sin 01 > |M|, and behaves as a normal insulator 
otherwise. 

The low-energy behavior of electrons in the Haldane 
model is governed by the dispersion relation near the two 
inequivalent Dirac points K and K' . We can expand the 
Hamiltonian to the linear order in the relative momentum 
k with respect to the two Dirac points, and obtain 



H (fc) = VF {TzkyCTx - k^ay) + y 



gTzCfz 



(4) 



Here, vf = 3at/2, A/2 = M, g = — 3\/3t2 sin and a 
constant energy — 3^2 cos has been omitted. CTq repre- 
sent two sublattices and = ±1 correspond to two Dirac 
valleys. Using a unitary transformation H' = U^HU 
in Eq. (4) with U — [1 4- CTj, + (1 — (Jy) Tz] /2, we obtain 
H' (fc) = Vf {kyCFx — k^cFy) -\- ^TzCFz -\- gcfz- To the linear 
order of fc, H'{k) is identical in form to Eq. (1), even 
though parameters in the two Hamiltonians have differ- 
ent physical meanings. More importantly, the condition 
for the occurrence of the quantum anomalous Hall effect 
in the Haldane model H, 3V3|i2sin(/.| > |M|, is also 
identical to that for a 3DTI thin film, i.e., \g\ > |A|/2 as 
shown in Fig. 1, due to the mapping of parameters below 
Eq. (4). Therefore, a 3DTI thin film in the presence of an 
exchange energy is a natural realization of the Haldane 
model. The exchange energy plays the role of the stag- 
gered fluxes in the Haldane model, and the hybridization 
gap of the surface states is equivalent to the alternating 
on-site energies. 

Quantum Hall Effect of 3DTI Thin Films. -We now 
study the QHE in a thin film of a nonmagnetic 3DTI. As 
a perpendicular magnetic field Bq is applied to the film, 
the action on the spin degrees of freedom is the Zeeman 
energy, which can be described by the last term, gaz in 
Eq. (1) (with g = g^f f ^.bBq)] and the action on the 
orbital motion can be included by use of the Peierls sub- 
stitution k — (k — eA) with A the vector potential due 
to -Bo- Since the unitary transformation U introduced 
above does not affect the orbital degrees of freedom, it is 
apparent that Hamiltonian (1) is still equivalent to the 
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Haldane model given by Eqs. (3) and (4) in this case, pro- 
vided that the same vector potential A is included into 
these equations. As a result, the QHE in a thin film of 
a 3DTI is equivalent to that in the Haldane model. The 
LLs in the Haldane model Eq. (4) subject to a magnetic 
field have been solved in the original work of Haldane [l| 



Er^.n = sgn(n) 



for nonzero integer n, and 

^,^,o = [.g+(A/2)T,]sgn(eBo) , 



(5) 



(6) 



for n = {]. Here, wi = ■\/2|e_Bo|f_F is the width of the 
V = 1 Hall plateau at = A = 0. When both g and A 
vanish, Eqs. (5) and (6) reduce to the standard LLs for 
massless Dirac fermions, which are additionally degener- 
ate for Tz — ±1, i.e., -E+,n = E^^n- Nonzero g and A 
may shift the relative positions of the LLs and cause the 
LLs to split, which will determine the quantization rule 
of the Hall conductivity. 

The Hall conductivity will be calculated numerically 
from Eq. (3) by setting — 7r/2 in a system of size 
64 X 64. The vector potential of the applied magnetic 
field is introduced into Eq. (3) via the integral form of 
the Peierls substitution t — te*^*^ (^2 — i2e*^'0j where 
9ij = e A-dl with the integral being along the electron 
hopping path. The magnetic fiux per hexagon is given by 
ip = % = 7r_Bo3-\/3a^/(/)o = 27r/Mo, where the sum- 
mation runs over six links around a hexagon, (po — h/e 
is the flux quantum, and Mq is an integer. The Hall con- 
ductivity at zero temperature can be calculated by use 
of the standard Kubo formula 



4e2 ^ 
S ^ 



em<Ep<e-a 



lu\{{n\v^ \m) {m\vy \n)) 



(7) 



Here, and e„ are the eigenenergies corresponding to 
occupied state |to) and empty state |n), respectively. S 
is the area of the hexagon lattice, and Va is the velocity 
operator with a = x oi y. With tuning chemical poten- 
tial Ep, a series of integer-quantized plateaus of Uxy will 
appear, each one corresponding to Ep moving in the gaps 
between two neighboring LLs. When g = A = 0, Eq. (3) 
reduces to the tight-binding model for graphcnc without 
spin degrees of freedom. It follows that the Hall conduc- 
tivity in the 3DTI thin film is half of that in graphene. 
The Hall conductivity for the 3DTI thin film is thus odd- 
integer quantized a^y — 2(^-1- with the prefactor 2 
originating from two surfaces, as expected. 

In Fig. 2, the calculated Hall conductivity a^y in units 
of ^ is plotted as a function of Ep/wi with magnetic 
flux ip = 27r/1024 for some different values of g and A. 
From Eqs. (5) and (6), we can see that for g ^ and 
A = 0, while their positions shift with g, all the LLs 
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FIG. 2: Calculated Hall conductivity Oxy in units of /h as 
a function of Ef with magnetic flux ip = 27r/1024 for (a) 
g = 0.5it)i and A = 0, (b) g = Q and A = wi, and (c) 
g — 0.5wi and A = 0.6iui. The corresponding LL energies 
for Tz = ±1 are shown in the right panel. 
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FIG. 3: Phase diagram for the QHE on the Ep/wi versus 
A/ Ac plane. The phase boundary is determined by the en- 
ergies of the LLs as functions of A/ Ac, calculated from Eqs. 
(5) and (6). Here, Ac — Wi/2g, and the Zeeman energy is 
taken to he g — O.Swi. 



are still degenerate for — ±1, so that the total Hall 
conductivity remains to be odd-integer quantized a^y = 
{2i + 1)^, as shown in Fig. 2(a). For g = and A ^ 0, 
the LLs with n ^ are degenerate for = ±1, but there 
is a splitting of the n — LL, yielding a new plateau of 
axy = 0, as shown in Fig. 2(b). For g ^ and A ^ 0, the 
additional degeneracy of all the LLs is generally lifted. 
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and each degenerate LL splits apart into two so that ah 

2 

integer Hah plateaus a^y = appear, as shown in Fig. 
2(c). 

In Fig. 3, we plot the energies of the LLs calculated 
from Eqs. (5) and (6) as functions of A for a fixed Zeeman 
energy, which effectively determine a phase diagram for 
the QHE on the Ep versus A plane. The splitting of each 
LL increases with A, which can be understood from Eqs. 
(5) and (6). Interestingly, when A is integer multiples 
of certain critical value Ac, the LLs for = ±1 and 
different n cross each other at the same time, leading to 
disappearance of nearly half of the Hall plateaus. For 



example, E. 



+ ,n 



E, 



,n+l 



at A = A,, E. 



E_ 



,n+2 



at 



A = 2 Ac, and on analogy of this, for n > 0. The LLs 
for n < cross in a similar manner. It is straightforward 
from Eq. (5) to obtain gAc = w\l2 or Ac = w\l2g. The 
Hall conductivity as a function of Ep is shown in Fig. 4(a) 
at A = Ac and 4(b) at A = 2Ac, respectively. In the 

2 

former, the Hall conductivity <Jxy = y\ is quantized into 
even-integer v = 2t plateaus plus a single odd-integer v ~ 
— 1 plateau, and in the latter, it displays odd-integer v — 
(2£-|-l) plateaus plus two even-integer and —2 plateaus. 
In both cases, the split Hall plateaus merge again partly, 
and most odd-integer [even-integer] plateaus disappear 
in Fig. 4(a) [4(b)]. 
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FIG. 4: Calculated Hall conductivity (Jxy in units of /h as 
a function of Ef with magnetic flux tp = 27r/1024 for (a) 
g — 0.5wi and A = wi, and (b) g = 0.5wi and A — 2wi. The 
corresponding LL energies for Tz — ±1 are shown in the right 
panel. 

Summary. - In summary, we have shown that the Hal- 
dane model for the QHE without LLs can be realized in 
condensed matter systems by use of a 3DTI thin film with 
an exchange field. As a perpendicular magnetic field is 



applied to the 3DTI thin film, Hamiltonian (1) and Eq. 
(4) for the Haldane model are still equivalent to each 
other, provided that the vector potential A is included 
in both equations. A rich phase diagram for the QHE in 
the 3DTI thin film is predicted as a consequence of the 
interplay between g and A. The simultaneous presence of 
nonzero g and A causes splitting of original odd-integer 
Hall plateaus axy = {2£ + l)x i^to axy = Remark- 
ably, when the product of g and A is at certain critical 
values, the plateaus can merge again partly, and most 
odd-integer (or even-integer) plateaus disappear. 
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